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Abstract 

Q_( A numerical study of the nonlinear wave solutions of the Rosenau-Pikovsky K{cos) 

equation is presented. This equation supports at least two kind of solitary waves with 
compact support: compactons of varying amplitude and speed, both bounded, and 
kovatons which have the maximum compacton amplitude, but arbitrary width. A 
new Pade numerical method is used to simulate the propagation and, with small ar- 
tificial viscosity added, the interaction between these kind of solitary waves. Several 
numerically induced phenomena that appear while propagating these compact trav- 
elling waves are discussed quantitatively, including self-similar forward and back- 
ward wavepackets. The collisions of compactons and kovatons show new phenomena 
such as the inversion of compactons and the generation of pairwise ripples decom- 
posing into small compacton-anticompacton pairs. 
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1 Introduction 



Compactons are compactly-supported solitary waves that appear as solutions 
of some generalizations of the Korteweg-de Vries (KdV) equation character- 
ized by the presence of nonlinear dispersion. Compactons were introduced for 
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the first time by Rosenau and Hyman [1] as solutions of the K(m, n) equation, 
given by 
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where x is the spatial variable and t is time, for certain values of m and n, 
and generalized by Khare and Cooper [2] for n = m. This kind of solitary 
waves presents some fundamental similarities with the solitons of the KdV 
equation, including the almost elastic interaction between them, but the non- 
linear dispersion in the K(m, n) equation makes compactons quite a distinct 
mathematical object with its own properties [3,4], being its compact support 
the most important. 

Numerical analysis has proved to be essential in the study of equations with 
compactly supported solitary waves. Pseudospectral [1,5], finite differences [6,7], 
Petrov-Garlerkin finite elements [8,9,10], Pade [11,12], and particle methods [5] 
have been used with more or less success. However, the numerical simulation 
of nonlinear travelling waves presents several numerically induced phenomena 
such as spurious radiation, artificial dissipation, or errors in group velocity, 
and the numerical evolution of compactons, especially at the edges of the sup- 
port, is not free of these undesired effects [8,13]. Special attention has to be 
paid to the behaviour of compactons under mutual collisions, a key property to 
check for their robustness [13,14,15]. Although there is a large number of non- 
linear partial differential equations presenting compactly-supported solitary 
wave solutions [16], only a few have been numerically studied in detail. 

Pikovsky and Rosenau [17] introduced a generalization of the K(n,n) equa- 
tion for the study of the quasicontinuous approximation of a one-dimensional 
lattice of coupled limit-cycle oscillators. Let us summarize the derivation pre- 
sented in Refs. [18,19]. The k-th self-sustained oscillator with frequency u is 
characterized by a phase <pk that obeys dipk/dr = u, where r is time. Under 
weakly coupling, the equations for the phase differences read 




= u + - (pk) + q(<Pk+i - <f>k), 



where q is the coupling function, even and 27r-periodic in its argument. In- 
troducing the variable u^ = 4>k+i — 4>ki the phase equations can be written 

as 
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whose quasicontinuous approximation yields 
or dv 3! dv 6 v ' 



where h is the spatial step of the lattice. Neglecting higher order terms and 
applying a change of variables, this equation can be non-dimensionalized to 

du dq{u) d 3 q(u) , . 

dt dx dx 3 



For q(u) = — cos(u) we obtain the so-called if (cos) equation, given by 

du dcos(u) d 3 cos(u) „ . 

1 — H r^- = 0, (4) 

dt dx dx 3 



that admits compactly-supported travelling waves such as compactons, kinks, 
and what Pikovsky and Rosenau referred to as kovatons, which result from 
gluing together a kink-antikink pair with a plateau of arbitrary width. 

The closed form expression for the shape of compactons and kovatons of the 
if (cos) equation is not known in closed form [18,20]. Its determination by 
means of numerical methods requires the solution of a nonlinear eigenvalue 
problem [21]. Recently, Garralon and Villatoro [22] proposed a procedure for 
the numerical quadrature of the solitary waves of the if (cos) equation, which 
can be used as initial conditions for compactons and kovatons numerical evolu- 
tion. These solutions are only characterized by one parameter, the velocity c of 
the solitary wave. The nonlinear relation between the amplitude and velocity 
of compactons must be numerically determined. 

The behaviour of compactons (and anticompactons) of the if (cos) equation 
differs from those of the if (2, 2) equation in several aspects [22]. Firstly, the 
absolute value of their velocity is upper bounded by c max = 2/%, the velocity 
of the kovatons (and antikovatons), whereas for the if (2, 2) there is no upper 
bound. Secondly, compactons (anticompactons) propagates to the left (right), 
i.e., to minus (plus) infinity in x, and the opposite occurs for those of the 
if (2, 2). And thirdly, the width of the if (cos) compactons increases as it does 
the absolute value of the velocity, as shown in Fig. 1, while it is independent 
of the velocity for the if (2, 2) ones. 

The if (cos) equation has several applications in physical problems such an 
array of superconducting Josephson junctions [23] and the spontaneous otoa- 
coustic emissions in the inner ear of certain lizards [24]. Let us also note that 
compact kinks are the solutions of nonlinear evolution equations modelling 
several applications such as waves in viscoelastic solids [25], double-stranded 
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Fig. 1. Compactons of the K(cos) equation with speeds c = 0.2, 0.4, 0.6, and 2/ir, 
and a kovaton with plateau length m = 4.0. 

DNA molecules [26], and nonlinear electrical transmission lines [27], to men- 
tion only a few. 

The contents of this paper are as follows. Next section presents a new numer- 
ical scheme used to solve the if (cos) equation, based on Pade approximants 
in space and a method of lines in time. Both compactons and kovatons have 
been studied, analyzing the residual structures that appear when propagating 
them alone as well as after pairwise interactions in Section 3. Finally, the last 
section is devoted to some conclusions. 



2 Numerical method 

Let us consider a method of lines in time and a Pade approximation method in 
space with periodic boundary conditions [12] for solving Eq. (4) numerically. 
Taking a fixed grid spacing Ax, the M+l spatial nodes of the interval 16 [0, L\ 
are Xj = j Ax, for j = 0, 1 . . . M, and the Pade method is written as 



where Uj(t) ~ u(xj,t), E is the shift operator, i.e., EUj = Uj+i, and the first 
and third derivatives are rationally approximated by means of B(E)/A(E) and 
C(E)/ A(E), respectively, where 
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(5) 



dt 



A(E) 



E" 2 + 26 E^ 1 + 66 + 26 E 1 + E 2 
120 



-E^-IOE^ + IOE^E 2 



C(E) 



24 Ax 

-E" 2 + 2E~ 1 -2E 1 + E 2 
2 Ax 3 ' 



The numerical method (5) is fourth-order accurate in space for regular enough 
solutions (u(x,t) G C 7 ), since its truncation error terms are given by 

TETMM)} = _|LJL «*(«) + 0(Ax 6 



It should be noted that the Pade numerical method (5) is a fourth-order 
approximation in Ax to the continuous equation (4), which in turn is a fourth- 
order approximation in h to the discrete equation (2). Although, Ax and h 
can be related by a proper scaling, physically it has no sense, since h is a 
fixed parameter of the original discrete model, but Ax is a tunable numerical 
parameter. 

In solutions of the iT(cos) equation with multiple colliding compactons and 
kovatons, ripples (or nonsmooth solutions) are developed reducing the effec- 
tive order of accuracy and introducing numerical instabilities which may blow 
up the solution [8,11,13]. In order to avoid these instabilities, an artificial 
viscosity term, <9 4 w/<9x 4 , with \x small enough, must be introduced into the 
non-dissipative method given by Eq. (5) and numerically discretized by means 
of a second-order accurate five-point difference formula, given by 

TKVMT E- 2 -4E- 1 + 6-4E 1 + E 2 rr 

V ( E ) U j = Ax 1 r ® 

Note that the addition of the artificial viscosity reduces the order of consis- 
tency of the method down to the second, resulting in 



TET„Hx,t)} = - r — +0(Ax 4 ). 



The second-order implicit midpoint rule is used for the integration in time, 
yielding 



jjm+l TTm / rjm+X , jjm \ 

A(E) i - i +(/iP(E) + c g(E)) ^ i 2 + i j 

/jjm+l + U m\ 

+ (B(E) + C(E)) cos 3 - = 0, (7) 
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Table 1 

Numerical errors in infinite norm for the one-compacton solution of the if (cos) 
equation at t = 20 as a function of Ax and fixed At = 0.001 (top), and as a 
function of At and fixed Ax = 0.1 (bottom), using L = 150, ti = 0, and cq = c. 

where U™ is an approximation to u(xj,t m ), with t m = mAt. 
3 Presentation of results 

Extensive numerical experiments with the above scheme has proved to be good 
in accuracy and conservation properties for both fi = and, when required, 
a small properly chosen /x. Let us summarize the results in the following sub- 
sections. 

3.1 One-compacton solution 

Let us consider the numerical simulation of the propagation of one compacton 
of the if (cos) equation calculated by method (7) without artificial viscosity 
(ji — 0), stopped in the integration interval by setting c = c. 

Table 1 shows the error in infinite norm between the exact one-compacton 
solution (determined as described in Ref. [22] with machine precision) and 
the numerical one at t — 20, i.e., maxj \Uj(20) — u(xj,20)\, as a function 
of Ax (top) and At (bottom). For Ax > 0.01 and At = 0.001, Table 1 
shows that the numerical error decreases as Ax does for compactons with 
velocities c = 0.5 and c = 2/n (the highest one). For Ax < 0.01 the error is 
higher for Ax < 0.001 the solution blows up even for a stopped compacton 
(indicated as — in the table); apparently this behaviour is surprising since the 
CFL (Courant-Friedrichs-Lewy) condition (Ax/ At > c — Co = 0) is always 
met. However, the CFL condition for the numerically induced wavepackets 
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Table 2 

Numerical error for the second invariant I2 for the one-compacton solution of the 
-RT(cos) at t = 20, as a function of Ax and fixed At = 0.001 (top), and as a function 
of At and fixed Ax = 0.1 (bottom), using L = 150, [i = 0, and cq = c. 

of radiation, to be studied later in this section, does not holds, hence an 
instability develops. On the other hand, Table 1 (bottom) shows that the 
error for both compactons is practically independent of At, and no problems of 
stability due to the CFL condition arise. Finally, let us note that the behaviour 
of the errors for anticompactons is the same as that for compactons with the 
same parameters (except that Co = — c must be used to stop them in the 
numerical frame). 

The numerical method used in this paper preserves exactly the first invariant 
for the .ff (cos) equation with fi = 0, i.e., 

00 

ii = J u(x, t) dx, 
—00 



since by summing in space it may be easily shown that J2m u m = J2m u m- 
However, second invariant 



sm(u(x, t)) dx 



is not exactly preserved. This invariant has been numerically calculated by 
means of the trapezoidal quadrature rule. 

Table 2 shows the error of the second invariant for one-compacton solution 
at t = 20 as a function of Ax (top) and At (bottom), i.e., | / 2 (20) — -^2 (0) | . 
As shown in this table, the error grows as Ax decreases for compactons with 
both velocities c = 0.5 and c = 2/tt, except in the case with the lowest Ax 
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Fig. 2. Forward (left plots) and backward (right plots) radiations generated during 
a compacton propagation with c = 0.5, Ax = 0.05, At = 0.01, and \x = 0, at time 
t = 250 (top plots) and t = 500 (bottom plots), highlighting their self-similarity. 

for c = 2/tt. By contrast, Table 2 (bottom) shows that the error of the second 
invariant is practically constant, nearly independent of At. 

The numerical evolution of compactons of the K(cos) equation shows the 
presence of two numerically-induced small amplitude wavepackets, one prop- 
agating to the left, the forward radiation, and the other one to the right, the 
backward radiation (recall that compactons travel to the left). Extensive nu- 
merical experiments show that the velocities of the wavefront for both the 
backward and forward radiation are independent of Ax and At parameters. 
On the other hand, the amplitude of radiations emitted for compactons in- 
creases as so do Ax and c. 



Another worth noting fact is that the radiations shape is self-similar, as in the 
Rosenau-Hyman K(2, 2) equation case reported in Ref. [13]. Both the forward 
and backward wavepackets present this self-similarity regardless of parameters 
Ax and At. Figure 2 clearly shows this property with plots of the radiations 
at two different times, stretching the horizontal axis in order to highlight the 
self-similarity. 

Similar results have been obtained for compactons with different amplitudes 
in long time simulations, enlarging the computational interval in order to 
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Table 3 

Front speeds of both the backward (c&) and forward (c/) radiation relative to that of 
a compacton with speeds c = 1/2 and c = 2/tt. The parameters used are Ax = 0.05, 
L = 300 and Co = c. 

avoid that the wavepackets overlap due to the periodicity of the boundary 
conditions. 

Table 3 shows the front speed of the forward (c/) and backward (c&) wavepack- 
ets of radiation relative to that of two compactons with speeds c = 0.5 and 
c = 2/tt. These speeds have been calculated by using linear regression of the 
position of the wavefronts in time, calculated by using a threshold equal to 
the half of the maximum amplitude of the radiation in the integration interval 
(the same procedure used in Ref. [13]). Table 3 shows that both the forward 
(cf) and backward (q,) front velocities are nearly constant (independent of 
both At and Ax) being Cf ~ 5 c and c& ~ — c. 

Due to the fact that the forward speed c/ is five times faster than that of a 
compacton, care must be taken when setting parameters Ax and At in order 
to meet the CFL condition, so that Ax /At > Cf ~ 5 c. This could be the 
reason of the blow up in the simulations of the last column of Tables 1 and 2 
(top). 

3.2 One-kovaton solution 

Let us consider the numerical simulation of the propagation of one kovaton 
of the -fT(cos) equation calculated by method (7) without artificial viscosity 
(/i = 0), stopped in the integration interval by setting c = c. 

Table 4 shows the error in infinite norm between the exact one-kovaton solution 
and the numerical one at t = 150 as a function of the plateau size m. The error 
is nearly constant for m < 4, but grows linearly for m > 4; in fact, a linear 
regression shows that, in such a case, the error is equal to (4.1 m — 0.49) x 10 -6 . 
Table 4 also shows that the error in the second invariant is very small and 
nearly independent of m. 

The linear behaviour of the error for a kovaton with m > 4 shown in Table 4 
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Table 4 

Numerical errors in infinite norm for the one-kovaton solution of the i^(cos) equation 
at t = 150 and for the second invariant as a function of the plateau size m, using 
Ax = 0.05, At = 0.01, L = 220, /i = 0, and c = c max . 
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Fig. 3. Evolution in time of the numerical error during one-kovaton propagation for 
t = 13 (top left plot), t = 47 (top right one), t = 71 (bottom left one), and t = 103 
(bottom right one), with m = 60, Ax = 0.05 and At = 0.01, cq = 0, and \i = 0. 
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is due to the appearance of a peak in its backward kink whose amplitude 
is dependent only on the value of m, as illustrated in Fig. 3. The evolution 
of the error presents three main features occurring simultaneously with the 
propagation of the kovaton, whose edges correspond to the vertical, broken- 
and-dotted lines in Fig. 3. First, a square pulse appears on the flat-top of the 
kovaton (Fig. 3, top left plot), whose right front propagates from the forward 
kink to the backward one with the kovaton's speed (Fig. 3, top right plot) until 
the kovaton has advanced a distance equal to m (Fig. 3, bottom left plot); here 
on, the shape of the square pulse remains constant without changes (Fig. 3, 
bottom right plot). Second, a downward peak appears at the position of the 
backward kink (Fig. 3, top left plot), whose amplitude grows until reaching a 
maximum value when the square pulse has passed through the flat-top (Fig. 3, 
bottom left plot); from this point on, the shape of the peak remains constant. 
And third, another square-like pulse, but with dispersive fronts at both edges, 
appears behind the kovaton (Fig. 3, top right plot) and propagates to the 
right, growing in width until it equals the length of the plateau, when it leaves 
the backward kink (Fig. 3, bottom right plot). The conservation of the first 
invariant of the if (cos) equation implies that the total area of the peak is 
equal to the sum of the areas of both square pulses. 

Kovaton propagation is also accompanied by two self-similar wavepackets of 
radiation whose shape is very similar to that of the one-compacton solution 
with maximum speed (cf. Fig. 2). The front speed of the forward (c/) and 
backward (q>) wavepackets of radiation are nearly constant, independent of 
both At, Ax, and the width m of the plateau, being approximately equal to 
Cf w 5 c and c& ~ — c. 



3.3 Compacton and kovaton collisions 

The numerical evolution of compactons and kovatons interactions is a more 
subtle problem than that of propagating them alone due to the large slopes 
developed, and a stable method for one-compacton or one-kovaton evolution 
may become unstable in mutual interaction. In fact, without artificial viscosity, 
/jL — 0, the numerical method blows up during any kind of interaction between 
them. On the other hand, artificial dissipation greater than zero generates tails 
and lags distorting the original signal. Numerical experimentation has proved 
a value of \i = 10~ 3 to be adequate. 

Figure 4 shows a sequence of plots with a collision between two compactons, 
one with speed c\ = 1/2 (the left one in the first snapshot), stopped in the 
computational frame by taking cq = c±, and the other one with the highest pos- 
sible speed, C2 = 2/n (the right one in the same snapshot). Both compactons 
collide elastically reemerging with the same shape after the interaction with a 
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Fig. 4. Collision between two compactons with velocities c\ = 1/2 (the left one in 
the first plot) and C2 = 2 /it (the right one), using /i = 10~ 3 , Co = c\, Ax = 0.05 
and At = 0.01. 

small shift in their positions, the largest compacton suffering a positive shift 
to the left, and the smallest one a negative one to the right (let's compare the 
position of the smallest one between the top and bottom plots in Fig. 4). In 
addition, a residual (not shown in the plot due to its small amplitude) appears 
at the position of the interaction. 

Figure 5 shows a interaction between a compacton with speed c = 0.9c max 
and a kovaton with plateau length m = 40. The compacton passes through 
the rising kink of the kovaton, then it reappears reversed (second snapshot 
in Fig. 5), travelling like an anticompacton along the plateau (fourth snap- 
shot), until reaching the falling kink of the kovaton (sixth snapshot), when 
it reemerges with its original shape (seventh snapshot). Of course, this phe- 
nomenon can not be clearly observed if the size of the plateau of the kovaton 
is less than the width of the compacton. Note that the width of the signal 
during the collision is the sum of both the width of the compacton and the 
kovaton before the interaction. Also note the phase shift in the positions of 
both solitary waves after the interaction. 

Compactons correspond to homoclinic orbits in the phase plane (u(-,t),u x (-,t)), 
at the critical point (0, 0), and kovatons to heteroclinic orbits between the crit- 
ical points (0, 0) and (it, 0); the evolution in time of these orbits is a useful tool 
to study their interactions. Figure 6 shows the collision between a compacton 
with speed c = 0.2 and a kovaton with plateau length m = 40. The interaction 
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Fig. 5. Collision between a compacton with speed c = 0.9c max and a kovaton with 
plateau m = 40, using fi = 10~ 3 , cq = c, Ax = 0.05 and At = 0.01. 

has actually three definite parts, the first being the entrance of the compacton 
into the kovaton, the second being the travelling of the compacton along the 
plateau of the kovaton, and the third being the exit of the compacton from 
the kovaton. In the top left plot of Fig. 6 the phase plane shows two ellipses, 
the largest one corresponding to the heteroclinic orbit of the kovaton and the 
smallest one to the homoclinic orbit of the compacton, both before the in- 
teraction begins. Next, the compacton enters the rising kink of the kovaton, 
corresponding to the left inner loop in the phase plane near the critical point 
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Fig. 6. Sequence of the phase plane of a collision between a compacton with speed 
c = 0.2 and a kovaton with plateau m = 40, with n = 10~ 3 , cq = 2/V, x £ [0, 200], 
Ax = 0.05, and At = 0.01. 

(0,0) (top center). When the compacton reaches the plateau of the kovaton, 
high slopes develop at the left side the plateau, causing the appearance of a 
sharp peak in the phase plane near (tt, 0) (top right). Then the compacton 
is almost completely inside the plateau of the kovaton, corresponding to the 
right inner loop in the phase plane near (ir, 0) (middle left). When the entrance 
of the compacton has finished, the compacton travels reversed all along the 
plateau of the kovaton, a remarkable situation that is plotted by two orbits 
in the phase plane, the smallest one corresponding now to an anticompacton 
on the plateau, and the largest one to the kovaton (middle center). When the 
compacton gets to the falling kink of the kovaton, it recovers its original shape 
and the phase plane plots the right inner loop at (tt, 0) first (middle right), the 
sharp peak at low amplitudes shortly afterwards (bottom left), the left inner 
loop at (0,0) (bottom center), and the final situation with the two original 
orbits (bottom right). 

Another fact worth mentioning is the appearance of a zero-mass small am- 
plitude ripple after every collision either between two compactons or between 
a compacton and a kovaton. Figure 7 zooms it for two compactons with ve- 
locities ci = 0.5 and c 2 = 2/V. The first plot shows the situation previous to 
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Fig. 7. Sequence of zooms of the evolution of the ripple developed after a collision 
between two compactons, one with speed c\ = 0.5 (the left one in the first zoom) 
and the other one with speed c 2 = 2/ir (the right one in the same zoom), using 
H = 10" 3 , c = 0.5, Ax = 0.05, and At = 0.01. 

the collision, with the faster (wider) compacton at the right. In the next plot 
we can see the birth of the ripple, only with its negative half part developed. 
The following plots shows the evolution in time of this ripple, travelling to 
the right in the computational frame at a speed equal to cq, i.e., it remains 
at rest at the point where the collision begins. This ripple is similar to those 
appearing in numerical simulations of K(2, 2) compactons, except for the fact 
that there are no sharp fronts inside (see [11]). 

The residual appearing after a compacton-kovaton collision consists of two sep- 
arated ripples, as depicted in Fig. 8 (top left plot), the first one corresponding 
to the collision of the compacton with the rising kink of the kovaton, labelled 
as a\ and a 2 in the plot, and the second one corresponding to the interaction 
of the now anticompacton at the plateau of the kovaton with its falling kink, 
labelled as b\ and b 2 . Thus, aq (g^) appears during the interaction between 
the negative (positive) slope part of the compacton and the rising kink of the 
kovaton. After the first ripple is generated, the compacton travels reversed 
along the plateau, and then b\ (b 2 ) appears during the interaction between 
the positive (negative) slope part of the now anticompacton and the falling 
kink of the kovaton. It should be noted that the distance between the two 
ripples equals the length of the plateau of the kovaton. In the following plots 
of the same figure, as the evolution progresses, the positive (negative) semi- 
ripples break down and decompose into trains of positive (negative) pulses of 
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Fig. 8. Evolution of the ripples generated after a collision between a compacton 
with speed c = 0.5 and a kovaton with plateau length m = 40, using \i = 10 -5 , 
Ax = 0.05, At = 0.01, and cq = 0. The snapshots correspond to time t = 525, 1110, 
2850, and 4830. 

decreasing amplitude, all of them with a width of 4 it, corresponding to small 
amplitude compactons (anticompactons). 



3.4 Compacton- anticompacton collisions 



Collisions between compactons and anticompactons, even adding the artifi- 
cial viscosity term, have proved to be unstable resulting in blow up, although 
extensive numerical simulations have shown that the behaviour of the numer- 
ical method for the propagation of anticompactons alone in terms of the error 
and the invariants is the same as that of compactons with the same speed, 
discussed in Section 3.1. 

Figure 9 shows the collision between a compacton with speed c\ = 0.5 and an 
anticompacton with c<i = —0.5. The top plots are two snapshots previous to 
the collision. In the middle ones, zoomed in, unexpected high frequency phe- 
nomena that begin to occur at the external edges can be observed. Conversely, 
the solution is apparently smooth at the central part. The amplitude of these 
high frequency phenomena increases as the blow up time approaches as shown 
in the bottom left zoomed out plot. Finally, the last bottom right plot shows 



1(3 




Fig. 9. Snapshots of a compacton-anticompacton collision, both with \c\ = 0.5, at 
times t = 16 (top left plot), 18 (top right), 20 (middle left), 21.8 (middle right), 22 
(bottom left), and 22.1 (bottom right), using Ax = 0.05, At = 0.01, Co = 0, and 
H = 10" 4 . 

the blow up. 

Figure 10 shows the phase plane plots corresponding to the same snapshots 
shown in Fig. 9. The left homoclinic orbit (negative amplitudes) of the top left 
plot in Fig. 10 corresponds to the anticompacton, and the right one (positive 
amplitudes) to the compacton. The top right plot in Fig. 10 show that as the 
travelling waves approach each other, the slope of the solution increases at the 
position where the compacton and anticompacton first met. In the middle letf 
plot of Fig. 10 a sharp peak in the bottom part of the melted orbits can also 
be observed, corresponding to the high frequency phenomena that begin to 
develop at the external edges of both travelling waves shown in Fig. 9 (middle 
left plot). These high frequency phenomena increase at the external edges 
causing the peak to grow in the phase plane as depicted in in Fig. 10 (middle 
right and bottom left plots). The final bottom right plot makes evident the 
blow up. 

The blow up during the compacton-anticompacton collision is due to the de- 
velopment of dispersive shock-like phenomena of increasing amplitude which 
appears in the regions of the solution with negative slope. Figure 11 (left plot) 
shows the solution for \x = 10~ 4 a few time steps before blow up; the use of 
a larger artificial viscosity (right plot) simply delays the blow up. The four 
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U(x,16) 



U(x,18) 




Fig. 10. Sequence of phase plane plots of a compacton-anticompacton collision both 
with |c| = 0.5 at times t = 16 (top left plot), 18 (top right), 20 (middle left), 
21.8 (middle right), 22 (bottom left), and 22.1 (bottom right), using Ax = 0.05, 
At = 0.01, c = 0, and \i = 10" 4 . 



|i = 0.0001 n = 0.001 




25 30 35 40 45 50 25 30 35 40 45 



Fig. 11. A compacton-anticompacton collision both with |c| = 0.5 at time t = 65.1 
with [i = 10~ 4 (left plot) and 10~ 3 (right one), using Ax = 0.05, At = 0.01, and 
c = 0. 



shocks shown in Fig. 11 present high frequency dispersion whose amplitude 
grows up as blow up time approaches. Our numerical simulations show that 
these phenomena is robust to changes in the parameters of the numerical 
method, Ax and At, indicating its possible origin in the analytical behaviour 
of the solutions of the K(cos) equation. 
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4 Conclusions 



The propagation of kovaton and compacton of the K(cos) equation has been 
studied by means of a new Pade numerical method. The good accuracy of 
the numerical method has been assessed by means of the error in the so- 
lution and the invariants for the propagation of compactons and kovatons. 
For compactons, small forward and backward wavepackets of radiation with 
a clear self-similar shape have been reported; the speed of the wave fronts of 
both radiations have been determined, being constant and independent of the 
grid spacing and time step of the numerical method. For kovatons, the main 
sources of the numerical error in the solution are the appearance of a small 
residual at the falling kink, whose area depends on the length of plateau, and 
the emission of a backward, square-like pulse whose length is also similar to 
that of the plateau. 

The study of the collisions between compactly supported, solitary waves of the 
if (cos) equation requires the addition of small artificial viscosity for stabil- 
ity, as in the case of the K(n, n) equation. For compactons, the collisions are 
very similar in both equations, showing a small amplitude residual decompos- 
ing in compacton-anticompacton pairs. The interaction between a compacton 
and a kovaton with a large enough plateau shows that the compacton passes 
through the rising kink, then it travels reversed like an anticompacton along 
the plateau, until reaching the falling kink of the kovaton, when it reemerges 
with its original shape. After the interaction two small residuals are gener- 
ated, which decomposes in trains of compacton-anticompacton pairs. Finally, 
the collisions between anticompactons and compactons result in the blow up of 
the solution due to the grow up in amplitude of dispersive shock-like structures 
in the solution. 

Present results show the great interest of the numerical study of the interac- 
tions of compactly supported solitary waves of generalizations of the K(n, n) 
equation, as those presented in Ref. [16], since only a few have been numeri- 
cally studied in detail. Furthermore, the use of other numerical methods, like 
finite volume methods for nonconvex conservation laws, is promising since the 
presence of shock-like phenomena can be considered a common feature of all 
the evolution equations with nonlinear dispersion. 
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